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Abstract 

The usual development of the continuous-time random walk (CTRW) proceeds by assuming that the present is 
one of the jumping times. Under this restrictive assumption integral equations for the propagator and mean escape 
times have been derived. We generalize these results to the case when the present is an arbitrary time by recourse to 
renewal theory. The case of Erlang distributed times is analyzed in detail. Several concrete examples are considered. 
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1. Introduction 

In this article we study transition probabilities and mean exit times of continuous-time random walks 
(CTRWs). By this we understand a random process X(t) whose evolution occurs purely via jumps of a 
random magnitude AX{t) that happen at random times t„, n € Z where the "waiting times" At„ = tn — tn-i 
are independent and identically distributed (iid). CTRWs have revealed as an interesting tool to model a 
large variety of physical phenomena that undergo sudden random changes. From a Mathematical perspective 
the consideration and interest in those processes can be traced back to the seminal work of Kolmogorov [1] 
and Feller [5] . Applications to describe changes of stock markets due to unexpected catastrophes were first 
noted in the seminal work of Merton [3] , where it is assumed that inter-catastrophe times are exponentially 
distributed independent of the magnitude of the catastrophe, i.e., that catastrophes are driven by a compound 
Poisson process (CPP). CTRWs generalize in an important way the latter processes allowing for general 
distribution of the waiting times. Correlation between waiting times and jumps is also permitted. 

In Statistical physics CTRWs became popular after the work of MontroU & Weiss [3], and have been 
used to describe physical phenomena ever since. To list a few examples we note applications to earthquake 
modelling (e.g., [516] ). rainfall description [7] or to transport in disordered media (e.g., [8]). More recently, 
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the use of CTRWs has been advocated to give a microscopic, tick-by-tick, description of financial markets: 



Unfortunately, the Markovian nature of CPP does not extend to the more general CTRW. In spite of 
this ominous situation, CTRWs satisfy a pseudo-Markovian property, namely that the knowledge of the 
past-prior to the last jump provides no further information to determine the future evolution than merely 
knowing the state of the system at such a jump time — cf. Eq. ([4]). This fact explains why scholars have 
usually focused in studying the statistics of CTRWs right after a jump occurs. In this regard, a linear 
integral equation for the mean escape times off a given interval has been derived |12I13| . Similarly, the basic 
probability for the process to be found in a certain region, given the position at a jump time, satisfies a 
certain integral equation first derived by Weiss |15|4j . 

However, this setting — wherein the present must be one of the jumping times — does not cover the most 
general situation. The relevance of this fact is further stressed by noting that in several physical problems 
one may not be able even to decide if such arrival has occurred. Hence, the issue of how to generalize the 
aforementioned framework to arbitrary present times arises naturally. In this sense some extensions have 
been considered in the past [15|16j . where it is assumed that the process starts at t = and behaves in a 
different way in the time period prior to the first jump. Nevertheless, to our knowledge, there is no robust 
and self-contained development on this topic in the literature. This paper tries to fill this gap and considers, 
in particular, the determination of the propagator. We also show how to obtain the mean exit time. We 
find that, unless jumping times are exponentially distributed, the results of |15I12| must be corrected. The 
implications of this fact to the calculation of the correlation function of the process are obvious. We consider 
here the simpler uncorrelated case: the more general case when jumps and sojourn times are correlated 
requires new ideas and will be the subject of a future publication. 

The structure of the paper is the following. In section [2] we obtain the unrestricted propagator of a CTRW 
and relate it with that corresponding to starting at jump times. This connection involves a certain object 
whose distribution may be found by recourse to classical renewal theorv [17|18| . General correlation functions 
follow immediately. Section [3] is devoted to solving the integral equation for the after-jump propagator by 
means of a joint Fourier-Laplace transform. In section |4] we extend the results of Masoliver et al. [12j and 
obtain the mean escape time off intervals .starting at an arbitrary instant. In section [SJ concrete, explicit 
formulae are given in the case that waiting times have Erlang distribution. Recently, Erlang times have been 
the subject of much interest in the context of information traffic and phone-calls waiting times |19I20| . In the 
context of transaction orders in financial markets the appearance of this distribution can also be expected 
since it takes, at least, two arrivals (buy and sell orders) for a transaction to be completed. For further 
applications to ruin problems and insurance see [21j . 

To help a reader not familiar with measure theory we assume that all distributions involved have a density. 
However, we find that all results extend to a general situation. 



2. Fundamentals of CTRWs 

Recall that any realization of a CTRW is given by a series of step functions in such a way that X{t) 
changes at random times . . . , t_i, to, ti, t2, • ■ • , in, • ■ • , while it remains fixed in place between successive 
steps — see Fig. [TJ The interval between these successive steps defines a sequence of independent identically 
distributed (iid) random variables At„ = t„ — tn-i, the waiting times. The (random) change, or jump, at tn 
in the process is given by AX„ = X{tn) — X{tn-i). We assume that Z„ = (AX„, At„) defines a sequence 
of iid two-dimensional random variables having joint and marginal densities p{x, t), h{x) and "(/'(r) given by 



see I9I10I11I12I13TT4] . 



p{x, T)dxdT = Pr (a; < AX„ < a; -I- dx, r < At„ < r + dr^ 



(1) 



h{x)dx = Pr {x < AXn <x + dx) , tp{T)dT = Pr (r < At„ < t + dr) . 



(2) 
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Fig. 1. A sample path of the process X(t). 
2.1. The propagator and related quantities 



One of the basic objectives within the theory of stochastic process is to obtain the future evolution given 
the actual state of the system; in this regard the main object is the conditional density of X{t) given the 
present position X{r) = Xr,r < t or, in more physical terms, the propagator p{x, t\xr, r).lft = r + T where 
r > the propagator reads 

p(x, r + T\xr, r)dx = Pr < X(r + t) < x + dx\X{r) = Xr) ■ (3) 

Note that r, t and every t„ are clock times; this is unlike r which gives the measure of a time interval. In 
order to have a frame of reference, we must specify some time origin. For convenience we take it to coincide 
with one of the jump times, to = say. In some problems tg may be identified as the starting point of the 
process X{t), but this is just one possible occurrence — see Fig.[TJ 

Now, recall that a generic CTRW is non-Markovian, and hence knowledge of the past gives additional 
information to that already provided by the state of the system at present. An exception to this occurs 
when the present happens to be one of the jump times: transition probabilities starting at a jumping time 
tn, show a pseudo-Markovian property: 

Pr [x < X{tn + t) <X + dx\X{tn) = Xn,X{tn^i) = X„_i, • . • ) = 

Pr {x < X{tn + T) <x + dx\X{tn) = Xn) = n(x, t„ + t|x„ , t„)da;, (4) 
and, whenever r < tn is 

Pr (x < X{tn +t) <X + dx\X{tn) = Xn,X{r) = Xr) = Uix, tn + T\Xn,tn)dx, (5) 

which is the backbone of these renewal processes. Intuitively, once a jump happens the system forgets 
all the past previous to the jump time starting anew. It must be stressed that Il(x,tn + T\xn,tn) in the 
expression above could be interpreted as the propagator starting from the jumping time i„ — i.e., roughly 
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n(a::, t„ + r|a;„, i„) is interpreted as p{x, i„ + r|a;„, i„) whereas when we write p{x, r + rjxr, r) the time r is 
arbitrary. Further, it can be proven that n(a;,t„ + T|a;„,t„) must satisfy the homogeneity condition 

Il{x, tn + T\Xn,tn) = ^{x ~ X„, t|0, 0) = Ti{x - X„, r). (6) 

These properties imply that n(x,t) solves the following integral equation (see Weiss [15]): 

n(2;,t) = (5(x)(l - *(t)) + / / p{x\t')\i{x - x\t ~ t')dx' dt' , (7) 



J -oo 



where '^{t) = Pr (At„ < t). Note that a class of CTRWs for which this equation can be solved in a closed 
form has been given in [iSj . We elaborate on this in section [3] below. 

However, although the assumption that the present is one of the jumping times may be convenient, there 
is, unfortunately, no convincing reason to that effect, neither from a mathematical nor a physical point 
of view. The evolution of the system is then described by the more general object p(x,r + T\xr,r) and 
the issue of determining it arises in a natural way. Here we address this problem, assuming for the sake 
of simplicity, that jump times and jump-magnitudes, At„ and AX„, arc mutually independent, i.e., that 
p{x,t) = ip(t)h(x). The general, correlated case will be the subject of a future publication. 

To this end set in Eq. jS]) p = p^''-' + p^^-* where p^"^ {B, r + T\xr,r) is the probability of a transition from 
Xr to B = [x, X + dx] in the time interval (r, r + r] occurring with no jumps. Obviously 

p'-°'^{B,r + T\xr,r)dx = Pr (^Xr+r e B, iV(,.^,,+^] = = x,.^ = 6{x - Xr)Pr (Ni^,.,r+T] = 0)d.T. 

Here N(^r,r+T] denotes the number of jumps of X on (r, r + t]. We next evaluate the probability that a such 
transition occurs with one or more jumps: 

p^^\B,r + T\r,x)dx = Pr (^Xr+r e B,iV(^^^+^] > l\Xr ^ (8) 

To this end wc introduce the "excess life" Tr as follows: r + is the time at which the first jump past r 
occurs. Note that is random, in contrast with r above, which is a number. Obviously, r + Tr = tn for 
some n (see Fig. [T] above). By conditioning respect to this object one has, by the total probability theorem 
and the memoryless property (j4]) 

p^'^\B,r + T\xr,r) = 

f f Pi (Xr+r e B, N^r^r+r] > 1 ^t„ = x' , t„ = t'\ Pr (Xt^ € dx' , tn € dt' Xr = Xr) ■ (9) 

Note that we use the convenient notation Pr [Xt^ e dx') = Pr (X(t„) e [x' , x' + dx']) and so forth. A further 
simplification arises from the independence assumptions: 

Pr (^Xr+r £ B, N^r.r+r] ^ 1 = x' , tn = t'^ = JI{x — x\ r + T — t')dx, 

Pr (^Xt^ e dx', tn e dt' Xr = Xr^ = h{x' - Xr)dx' Pr (rr e [t' - r,t' - r + dt']^ . 
Hence, by substitution into the RHS of ([3]) wc find 

p{x, r + T\xr, r) = 6{x — Xr) Pr (r^ J J — x' ,t — t') Pr ^r^ S dr'^ h{x' — Xr)dx' = 

<5(a;-a;^)(l- $(r|r)) + / / U{x - Xr - x' ,t - T')(j}{T'\r)h{x')dx'dT' . (10) 

Jo J ~oo 

Here $(r|r) and 4>{T\r) denote accumulated distribution function and the density of r^: 

$(r|r) = Pr (t^ < t) ^ Pr (t„_i < r < i„ < r + t) and (/)(T|r) = Mllll. (n) 

n=l 

To summarize, the propagator is recovered by (fTO|) if $(T|r) and Ii{x,t) are known PI The latter can be 



^ Two-point and correlation functions follow immediately from Eq. IIIUI I. 
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completely retrieved by means of a joint Fourier-Laplace transform (see section [3]). It turns out that the 
former can also be recovered in closed form by means of classical renewal theory as we now show. To this 
end let m{t) denote the renewal function, the mean number of jumps between to = and t. It satisfies the 
integral renewal equation (see [17] or [18]) m{t) = '^{t) + 4'(t — t')dm{t'). Further, by the law of total 
probability one finds that $(T|r) solves the following renewal equation 

$(T|r) = *(r + t) +y (^$(r|r - r') - l)'i*(»''): (12) 

whereupon, appealing to the renewal's theorem one finds that ^{T\r), (j)(T\r) are given by 

$(r|r) ^ ^{r + t) - (^1 - *(r + t - t')^dm{t'), (l>{T\r) ^ ilj{r + t) + ip{r + t - t')dm(t'). (13) 

Laplace transformation is useful to evaluate in an explicit form the latter objects. To this end note that 
manipulation of Eq. (fT3|) yields also 

$(T|r)=^ (l-'^{r + T -t')^dm{t'). (14) 

Let 4>{s\r) and be the Laplace transforms of 4'{T\r) and ipi^)- Upon Laplace transformation, we find 

^'(t).^^^^r°°e-^d.,c>0, (15) 

1 i-c+ico I-CO 

0(r|r) = / e'^^+'Ul-iis)) e-'''dm{t)ds, c> 0. (16) 

StI"* Jc-ioo ^ ^ Jr 

The above results recover the propagator in closed form. We next mention properties of this object that 
follow from Eq. (jlOp . Note first that p is spatially-homogeneous: p{x,r + T\xr,r) = p{x — Xj.,r + T|0,r) = 
p{x — Xr, r + rjr). This in turn implies that Xf^ — X^^ is independent of Xt-^ for any two times ti < ^2- This 
does not imply, though, the stronger property of independence of the increments as in Levy processes. 

By contrast, p is homogeneous in time only when $(T|r) is independent of r: <I>(T|r) = $(r|0) = ^'(t). 
Similarly, a comparison between ^ and (fTO|) shows that for the condition p{x,r + rjr) = II{x,t) to hold 
one needs, again, that $(T|r) = <i>(T|0). It can be proven that the latter holds only if waiting times are 
exponentially distributed. If this is not the case the after-jump propagator 11 does not describe the temporal 
evolution of the system; instead, the more general object p{x — Xr, r + rjr) must be used. 

The non-Markovian character of the process (see section [3|) implies that any available information about 
prior evolution of the process affects the transition probability. Suppose for instance we wish to determine 
the future evolution of the system if we know that X{r) = x^ and, in addition, the value of the previous 
jumping time, say t„_i. With no loss of generality set t„_i = to = 0. The transition probability adapted to 
that information is 

Tr{x, r + T\xr, r)dx = Pr [x < X{r + t) < x + dx\X{r') = Xr^r' < r) . (17) 
With similar ideas to those used before one can prove that 

7r(.T, r + T\xr,r) = Tr{x — Xr, r + t|0, r) = 

5{X - Xr)^—^^^^^^ + r r h{x') ^^'' ^ ^'] lV{x ~Xr~x',T~ T')dx' dr' . (18) 

l-*(r) 7o J-oo 'l-l'(r) 

3. Transition probabilities in frequency domain 

Since Eq. ([7]) is of convolution type an explicit analytical solution can be retrieved in the Fourier-Laplace 
domain. We remind how this is done — for other considerations on this regard see also |15I13I22I23| . By taking 



5 



a Fourier (Laplace) transform in space (respectively, in time) Eq. ([7]) yields that the joint Fourier-Laplace 
transform of 11 is recovered as: 



OO /'OO 



n(a;,i)e*' 



dxdt = 



l-^(s) 



^1 - h{uj)i;{s))s 

Here w is a real variable while s = + isi is complex, > and we define 



(19) 



(20) 



By inversion one could, in principle, recover both p{x, r + rjr) and n(a;, r). Regarding the latter we note 
that n must be a combination of an isolated mass at x = and a sub-stochastic density n(^^(a:, s) which is 
the inverse FT of Q = n(cj, s) — n(oo, s). Thus, the inversion formula reads 



n(x,t) = 5(x)(i-*(t)) + ^ 



+ CXD 



(21) 



Recall that for Xt to be Markovian, both the propagator p and after-jump propagator 11 must satisfy the 
Chapman-Kolmogorov equation. Using the spatial invariance, we find that this equation reads, in Fourier 
domain, as n(a;, t + I) = II{lu, t)Il{u!, I) for all uj and times t, I. This is Cauchy's functional whose solution is 
Il{u!,t) ~ A(ijj)e^^*. Taking a further Laplace transformation (in time) a comparison with Eq. p9|) shows 
that we must require A{u!) = l,B{u!) = \(l — h{Lu)) and V'(s) = A/(A + s) for some positive A, i.e. that 
times must be exponentially distributed for the process to be Markov. 

For illustrative purposes, we consider an example. Assume that the jump probability is exponentially 
distributed with parameter 7 but with different probabilities to jump left and right: 



Kx) = 7(^ + i^)e ''^1:e>o + - K)e''^l^<o, 
where — 5 < k < h and 7 > are parameters. In this case we find n(i)(x-,s) = 



"f1p< 



1 - 2K'^tp 



1 - 2kV 



(22) 



(23) 



where = V'('S) ^nd (p = ^{s) = 2y 1 — V + n^ip'^- Inversion of the corresponding Laplace transform can be 
accomplished when n ~ ^, h{x) ~ 76 '^l^>o and 'ip{s) = jr^- In this case one finds that 



p{x, r + rjr) = n(a;, r) 



5{x) 



^e-^--^-/i(2v/A^)l,>o, 



(24) 



where /«(■) is the modified Bessel function of the first kind with order n. For general values of k one 
must resort to numerical techniques. In Fig. [2] we perform a numerical inversion of Eq. (|23p corresponding to 
= aTs different values of the parameter k. In the upper panel we plot the "density contribution" 

to the propagator n(^)(a;,r) = n(a;,T) — e~^'^5{x) as a function of 7a; for fixed value of time. Notice how 
the discontinuity of h{x) at a; = — cf. Eq. ((22|) — is inherited by Ii^^\x,T). The lower panel plots the 
accumulated distribution function F{x,t) = n(a;', t)c?x', whose discontinuity comes from the delta 
contribution to n(a;, r) at a; = 0. 

Returning to the general case, we now derive the propagator. Operating with the Fourier-Laplace transform 
and using (fT5|) . pB]) and we obtain that in the Fourier-Laplace domain, Eq. PT7)1 reads 

p(w, s\r) = (1 - 0(s|r)) Is + ^{s\r)h{Lo)i\{uj , s) (25) 

The representation shows that the transition probability contains an isolated mass at a; = — which gives 
the probability that no jump has taken place — and a sub-stochastic density p^^^ {x, r + rjr). Upon inversion 



p{x, r + T\r) = 5{x){l - <f>(t|r)^ -^p'-^^x, r + r|r). 




Fig. 2. Numerical inversion of the "density contribution" to the propagator (upper panel) and the accumulated distribution 
function (lower panel) for exponential waiting times with intensity A. The parameters were chosen such that Ar = 2.0. We 
consider several possible values of the bias-related parameter k, from the symmetric case fc = to the pure persistent case 



1 /'C+ioo p-\-oo 

where p^^^x, r + r|r) = — / / e"*-*"^0(s|r)/i(w)II(a', s)dujds, c> 0. (26) 

47r I J c—ioo J — oo 

Eq. ([3]) gives the propagator in a fully explicit way. We skip a similar expression for Tr(x,T + r\r). 

The above expressions simplify if r — > oo, which corresponds to the steady-state solution. This limit is 
relevant both from the mathematical and the applied point of view, since it describes the case in which the 
only information available to the observer is the present value of the stochastic process. Using equation (|25[) 
and setting 4>oo{t) = lim 0(r|r) the full propagator reads 

r — ^oo 

lim p{u^, s\r) = L^22^^ + V / ^ (27) 

Eq. ((TT]) coincides formally with an equation already reported in [TB] . We admit that the resemblance is not 
spurious, but it must be stressed that the problem addressed there, although related, is different from ours. 
Nevertheless, we can go a step further by evaluating (jjooiT) explicitly; indeed, the renewal theorem yields 
that lim m(<) = t/ii, where M = /q ttp{t)dt is the mean sojourn time. Inserting this into Eq. we obtain 



0oo(t) = 



1 - 4'(r) 



1 /■°° 

$oo(t) = lim $(T|r) = 1 - - / td-^it + t). (28) 



4. Mean exit times 



Another significant magnitude in CTRW problems is the mean exit time T{xr,r) off a given interval 
A = [a,b]. This is defined as follows: suppose we observe that at a certain moment r the process was in 
Xr ■ X(r) = Xr, Xr G [a, b]. Then T[a^b]{xr, f) is the mean time that takes for the process to exit A. Note that 
for ease of notation we usually write T^a,b]{xr,'r) = T{xr,r). When r is one of the jump times, r = tn say, 
we can simply write T(a;„) = T[a,b]{xmtn). It can be proven that T(-) satisfies the following linear integral 
equation (see [12]) 

T{xn) = l-i+f h{x' - Xn)T{x)dx'. (29) 

J a 

We next show how to generalize these results to obtain T[a^b]{xr,r). Let r + r,. and r + t^'' be the first 
jump time and, respectively, the first time past r at which X{t) exits A so that T[a^b]ixnr) tf^ > — see 
section [2] and Fig.IH Obviously t^"^ — Tr if the process exits A at the first jump after r while tj;'' = t,. 

if the process jumps to a position y = Xr+Tr •''tm within A. Thus we can write i^'' = r,. + t'^^'l^^'' lx{r+Tr)eA- 
By taking an average we find that 

< tf;+;'- lx(r+r^)eA >= / < <r+r„ > (^r+r,. S dy\Xr = X, ) , (30) 
J a 

where, in the spirit of the total probability theorem, we have conditioned on the after-jump position X{r+Tr) 
and used the "lack of memory past jump-times" property, viz. Eq. ([4]). This expression is simplified further 
by noting that < t^.^^ >= T{y). Further, 

Pr {Xr+T^ G [y,y + dy]\Xr = a;^) = Pr {AX e [y - Xr.y - Xr + dy]) = h{y - Xr)dy, (31) 

since the jump amount is statistically independent of the previous position. Hence, the RHS of (|30p equals 

T[y)h[y — Xr)dy. By substitution we finally have 

T{Xr,r)=<Tr> -^L + T{Xr) = T{Xr)+ t' {^(l){t' \r) - iP [t' )^ dt' . (32) 

Thus the mean time to exit [a, b] past r is obtained by subtracting /.i— < > to the mean exit time 
after a jump happened, Tixr) — which solves the linear integral equation (|29[) . In section 5 (see Eq. (|4ip ) 
we evaluate this correction in the case of Erlang waiting times. 



8 



5. Erlang sojourn times 



5.1. Distribution and mean of the first jump past r 



Even although for a large variety of physical situations waiting times can be fitted by an exponential 
distribution, there exist other possibilities of interest. Here we consider an alternative distribution. Recall 
that a random variable is said to have Erlang distribution £r{i', A) whenever it can be expressed as a sum of 
i^, G N, iid exponential variables with parameter A. It turns out that £r{i', A) coincides with the Gamma 
distribution r(/3, A) whenever /? = is an integer. This corresponds to 

.-At . ^ 



/-I e 



A + s 



(33) 



To recover 0(T|r) we first consider the case of the renewal function, cf. Eq. ^T5\i . Here the integrand has poles 
at points bj = X{ej — 1) where ej = ,j = 1, . . . , i/. Hence, by closing the contour appropriately in the 
complex s-plane and integrating by residues we find that the renewal function is (note that bjn = Re bj < 0) 



,-A(l-e,)t 



(34) 



Inserting this into (|16[) we find upon Laplace inversion that 



c^{r\r) = ±^l _ dse-Y. 



2m 



A6, 



1 - 



A + s 



Yl •^(^■^)" 7 TT^(^n{r) where a„(r) - V e"e''^''. 



Further, the identity X]j=i ~ ^^nv shows that Q!„(0) = and (/'(t|0) = ')/'(''")• Also 



<^(T\r) = 1 - ^ (At 



ri-l e 



-At I' 



(n-l)! ^ 

We next evaluate the mean of this distribution. The upper expression in Eq. (jS.ip yields 

d4)[s\r) 



t'(l){t'\r)dt' = -- 



In particular, if = 2, then 



<i>(T|r) 



-At 



-2Ar 



I' + 1 



,6, r 



/A. 



Ar 



< > = 



t'(f>{t'\r)dt' 



3 + e- 



-2Ar 



2A 



(35) 



(36) 



(37) 



(38) 



Some of these calculations can be extended to cover the more general Gamma distribution T{i>, A), whose 
density is still given by Eq. (j33p . whenever i/ is a rational number: v = ^ with n,q irreducible integers. 
Closing the contour appropriately in the complex s-plane we find the derivative of the renewal function as 

'^^j Jit , ^ '"° 



.'(t)=V^eM + _l, 



i=i 



d^eAt(.-i) I' 1 _ 



In some cases the above integral can be evaluated in closed form. For instance, if ;/ = 1/2 we will have 

m\t) = \/^e-^* + AErfc( - a/aI), m(t) = i/^e"^* + ^^i±lErfc( - ^/A^) - 1, 
\ nt V TT 2 2 

where Erfc(-) is the complementary error function. 



(39) 



(40) 



5.2. Propagator and escape times under the Erlang distribution 



Given that Xr = x the expected remaining time to exit [a, 6], a < x < b, foUows by inserting (|37|) into 
32)) : since the mean of Erlang distribution is fi = j we find that 



T(x.,r)=T(x.)+(i^ + Eu^ 



(41) 



In particular, if i/ = 2, then T[a,b]ixr,r) ~ T{xr) — ^1 



-2Xr 



/2A. 



We next consider the after-jump transition probability. Note that for Erlang times £r{i', A) Eq. (fTO|) yields 



11(10, t) 



2m 



c-\-ioQ 



ds 



(A- 



(42) 



The integrand has poles at sj ~ X{q{uj)€j — 1), j = 1, . . . , where g(w) is a determination of h}/^{Lo). Note 
that < < 1 and all poles are located on the left half-plane. Thus, Cauchy's theorem yields that the 
"after-jump" propagator is n(x, t) = 5{x) I 1 — 5'(<) ) + II'^' (cc, t), where 



E 



1 



(At; 



ri-l ■ 



e«g - 1 



(n-1)! 



When V = 1, ([^S]) is the well-known inversion formula for a CPP. When v = 2 one has that 



n(i)(2;,o 



27r 



-'ioja:; — At 



--i= sinh (^A VXi^ + cosh (^AVX<^ - 1 - Ai 



duj. 



(43) 



(44) 



Note that, 11'^^' (x, i) = O(t^) as t ^ 0, which confirms again that the process can not be Markov. 

In Fig. [3] we perform a numerical inversion of Eq. (|43p corresponding to the jump density h{x) of Eq. ([22]) 
where k, = 0.1 and we take different values of the parameter v, concretely v = 1,2,3,4. In order to keep 
the plots commensurable wc have changed the A parameter in such a way that the no-jump probability 
1 — 5'(t) remains fixed. Observe the different decay behavior which is noticeable both in n^^^(a;,T) and in 
the accumulated distribution function F(x,t). 

The propagator p(x, r 4- r\r) can be evaluated in closed form when r — > oo. For v = 2, say, upon inversion 
of ([?7|) with (32), we find that Poo(a^, t) = \im p{x,T + r\r) reads 



Poo(x,t) = (l - '^oc:{t)^S{x) + 



1 

2^ 



^—iujx — XT 



2\fh 



sinh ( A 



cosh ( AV^r^ — 1 — -^r 



dijj. (45) 



Notice how it docs not quite settle to the after-jump propagator n(a;, t) of Eq. (|44p . In particular for r ^ 0, 
Poo(x,t) exhibits a quite different behavior, since Eq. (|45p implies that 



Voo{x,t) - ( 1 - - 



h{x)_ 



(46) 



where we recall that /i = 2/ A. Wc point out that the expansion pS)) also holds in the general case as can be 
easily seen using Eq. ([77)) . 
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